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Abstract 

For  the  dissemination  of precise  time  and  frequency ,  the  use  of  the  Two-Way  Satellite  Time 
and  Frequency  Transfer  (TWSTFT)  method  has ,  in  recent  years,  become  increasingly  valuable . 
By  the  application  of  spread-spectrum  technology ,  a  client  station ,  located  anywhere  within  a 
common  satellite’s  footprint ,  can  link  to  a  reference  station  and  compare  or  synchronize  its 
clock  to  the  reference  clock  within  nanoseconds  or  better  with  high  levels  of  confidence .  But  as 
the  frequency  stabilities  of  today’s  sophisticated  atomic-level  clocks  improve,  so  must  the 
stability  of  a  Two-Way  system’s  method  of  time  and  frequency  transfer.  Carrier-phase 
information  holds  the  promise  of  improving  the  stabilities  of  TWSTFT  measurements,  because 
of  the  great  precision  at  which  frequency  transfers  can  be  achieved.  The  technique  requires 
that  each  site  observe  both  its  own  satellite-translated  signal  and  that  of  the  cooperating  site.  As 
was  reported  in  earlier  papers,  the  translation  frequency  of  the  satellite  itself  is  an  additional 
unknown  factor  in  the  measurement,  and  it  must  be  taken  into  consideration.  However,  it  can 
be  shown  that  simple  estimates  of  the  satellite’s  local  oscillator  (LO)  frequency  will  suffice. 
Recent  work  has  been  conducted  with  implementation  of  the  signal  carrier-phase  in  operational 
TWSTFT  links.  The  purpose  of  this  paper  is  to  discuss  this  recent  work,  with  emphasis  on  the 
system ’s  development. 


I.  INTRODUCTION 

At  the  present  state  of  the  art,  one  of  the  primary  tools  for  the  dissemination  of  precise  time  and  frequency 
is  a  Two-Way  Satellite  Time  and  Frequency  Transfer  (TWSTFT)  system.  With  such  a  system,  a  client 
station  and  a  reference  station  each  employ  a  code  division  multiple  access  (CDMA)  scheme  to  lock  onto 
each  other  and  establish  communication.  The  spreading  function  used  for  the  CDMA,  while  providing 
the  spread-spectrum  capability  and  all  of  the  associated  benefits,  also  provides  time  markers  from  each 
station’s  respective  timescale.  This  allows  the  client  station  to  synchronize  its  time  to  that  of  the  linked 
reference  station.  Such  a  system  is  illustrated  below  in  Figure  1  [1]. 
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Figure  1.  A  basic  Two-Way  system. 


A  CDMA-based  TWSTFT  system  measures  the  time  difference  between  two  clocks,  and  can  provide 
frequency  by  differentiating  the  time-difference  measurements.  However,  the  CDMA-based  TWSTFT 
method  utilizes  the  carrier  wave  of  the  signal  and  not  the  superimposed  lower-frequency  pseudorandom 
noise  (PRN)  code.  Therefore,  CDMA-based  TWSTFT  suffers  more  from  multipath  effects,  which  is  a 
function  of  effective  wavelength,  than  does  the  carrier  signal.  Additionally,  it  is  only  the  change  of  the 
relative  phase  that  matters;  an  initial  phase  offset  induced  by  multipath  and  the  ionosphere  has  no 
relevance.  For  the  same  reason,  other  time  delays  introduced  into  the  system,  such  as  Sagnac  effects  and 
uncalibrated  static  time  delays  within  user  equipment,  will  have  less  effect  on  a  carrier-phase  system  than 
in  the  code-based  TWSTFT  configuration. 

Previous  efforts  to  use  the  carrier-phase  in  TWSTFT  have  developed  the  system’s  fundamental  equations, 
and  the  results  have  shown  some  promise  [2].  In  this  paper,  we  discuss  the  derivation  of  the  system’s 
fundamental  equations,  followed  by  the  results  of  two  test  runs:  a  collocated  test  at  USNO,  and  an 
intercontinental  test  between  PTB  and  USNO. 


II.  CREATION  OF  THE  SOLUTION  EQUATION  SET 

In  general,  the  carrier  frequency  of  a  received  signal  will  be  offset  from  the  frequency  generated  by  the 
transmitter.  The  magnitude  of  the  frequency  offset  can  be  mostly  attributed  to  four  system-variables 
(illustrated  below  in  Figure  2):  fSLO  ,  ,  k2 ,  and  df  [2],  which  are  defined  in  the  discussion  to  follow. 
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Figure  2.  Carrier-Phase  Two-Way. 
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As  shown  in  the  above  figure,  the  reference  station’s  (Station  1)  time  is  based  on  a  frequency  fsys  and  the 
frequency  of  Clock  2  is  offset  from  fsys  by  a  value  df  .  Realization  of  the  offset  between  the  frequencies 
of  Clock  1  and  Clock  2  is  the  basis  of  this  paper. 

A  simple  comparison  of  the  transferred  signals  will  not  yield  a  true  value  of  df  due  to  three  frequency 
perturbations  (or  frequency  shifts)  to  which  the  carriers  are  subjected:  two  caused  by  the  motion  of  the 
satellite  and  one  caused  by  the  satellite’s  unknown  local  oscillator  frequency,  or  fsw.  As  an  example, 
consider  a  one-sided  carrier-frequency  transfer  depicted  in  Figure  2,  where  Station  1  is  transmitting  a 
signal  to  Station  2.  The  signal,  which  is  relayed  to  Station  2  by  a  satellite,  is  subject  to  the  Doppler 
effects  resulting  from  the  slight  motion  of  the  satellite.  The  first-order  Doppler  coefficients  are  described 
by  the  following  two  equations. 


v 

*1=^  (1) 

c 

k2=I±-  (2) 

c 

where  Vn  is  a  projection  of  the  velocity  vector  of  the  satellite  in  the  direction  of  Station  n  and  c  is  the 
speed  of  light  through  the  transmission  medium.  Note  that  in  Figure  2,  the  velocity  of  the  satellite  has 
been  represented  as  a  vector  Vs  .  Suppose  that  Station  1  transmits  a  signal  with  a  frequency  flx ,  where 
fu  is  defined  with  respect  to  (w.r.t.)  the  system  frequency  fsys .  At  the  satellite,  the  frequency  of  the 
received  signal  is  centered  at  the  Doppler-shifted  frequency  ftx  +  k]ftx ,  which  is  then  mixed  with  the 
satellite’s  local  oscillator  (LO)  and  retransmitted  at  the  carrier  frequency  given  by  the  following  equation: 

fsv  ~  ftx  +  Kftx  —  f SLO  •  (3) 

The  carrier  frequency  is  then  observed  Doppler-shifted  at  the  downlink  site  as: 

fd  =  (/te  +kJtx-fsw)  +  k2(f,x+kJtX  ~  fsLO )•  (4) 

All  of  the  frequencies  given  in  Equation  (4)  are  referenced  to  the  system  frequency.  However,  the  value 
of  fd  will  be  measured  with  reference  to  the  receiver’s  local  clock.  Equation  (4)  must,  therefore,  be  re¬ 
written  such  that  its  terms  are  expressed  with  reference  to  the  clock  local  to  Station  2. 


Figure  3.  System  frames  of  reference. 
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For  the  purpose  of  converting  a  given  frequency,  as  defined  by  a  frame  of  reference  m,  to  its  equivalent 
frequency,  in  a  frame  of  reference  n ,  a  transformation  equation  has  been  derived: 


/„=/' 


V  Mn 


(5) 


where  fm  is  the  value  of  the  frequency  (w.r.t.  the  transmitter’s  clock),  which  is  to  be  transformed;  fn  is 
the  value  to  be  determined;  frefm  is  the  instantaneous  frequency  of  the  reference  clock  at  the  transmitter 
(given  w.r.t.  system  time);  and  frefn  is  the  instantaneous  frequency  of  the  reference  clock  at  the  receiver 
(given  w.r.t.  system  time).  Applying  this  transformation  equation  to  the  frequency  transfer  given  by 
Equation  (4)  produces 
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where  the  notation  Fmn  is  used  to  denote  a  frequency  that  has  been  transmitted  by  Station  m  and  measured 
by  Station  n.  Similarly,  we  may  derive  the  frequency  transfer  from  Station  2  to  Station  1  as  follows: 


^21=/, 


\  f sy*  J 


(l  +  *1)(l  +  *2)-/ffl0(l  +  *1). 


(V) 


In  addition  to  the  two  transfers  given  by  Equations  (6)  and  (7),  each  station  may  receive  its  own  signal. 
These  are  modeled  by  the  following  two  loopback  equations: 

Fn=ftx(^kl)2-fsLo(^kl)  (8) 


^22  —  fix  (1  +  ^2)  f  SI 


fsv 


Kfsys  +df  j 


(1  +  ^2  ) 


(9) 


With  theoretically  ideal  measurements,  the  values  Fn ,  F22 ,  F]2 ,  and  F2l  are  independent;  however,  under 
the  limitation  of  the  system  equipment’s  quantization  noise,  the  following  equality  is  satisfied: 

P\\  +  P22  ~  F\2  ^21  (10) 


This  equality  allows  for  an  immediate  data  integrity  check.  Alternatively,  it  will  be  shown  that  it  also 
allows  the  user  to  solve  the  system  using  only  three  of  the  four  equations,  thus  sparing  a  receiver  channel. 

With  the  two  carrier  transfers  and  the  two  loopbacks  defined  above,  the  carrier-phase  system  is  fully 
described  and  we  now  have  four  equations  from  which  we  may  resolve  our  four  system  variables. 
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III.  UNDERSTANDING  THE  SYSTEM  EQUATION  SET 

Generally  speaking,  the  frequencies  Fu ,  F22,  Fn,  and  F2X  are  provided  by  the  two  spread- spectrum 
modems  employed  by  the  system,  and  the  frequencies  are  reported  in  the  IF -band.  Because  the  carrier- 
phase  system  is  based  on  these  measurements,  it  may  be  insightful  to  examine  the  construction  of  the 
values.  Consider  the  following  system,  where  the  satellite  is  stationary  with  respect  to  both  TWSTT 
stations: 
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Figure  4.  RF  transceiver  signal  propagations. 


Relating  the  above  figure  to  the  carrier-phase  system,  the  illustrated  signals  have  the  following  associated 
frequencies: 

CLK]  — »  fsys 
CLK2  — >  fsys  +  df 
TX  i— >/* 

TX2^fa+y0flx 

LOA  f LOtx 

LOb  fiorx 
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LOd  fwa  +  y  of  LOtx  ’ 


where  the  instantaneous  fractional  frequency  y0  is  defined  as 


y  o  = 


dj_ 

f 


(11) 


As  shown,  the  RF/IF  mixers  are  supplied  tunable  synthesized  frequencies  from  their  respective  LO’s.  It 
can  be  shown  that  the  output  frequencies  of  the  synthesizers  will  have  the  same  fractional  frequencies  as 
their  inputs,  with  regard  to  fsys .  Now  if  the  signal  is  traced  from  left  to  right  (i.e.  from  the  signal  in  the 

form  of  TXx  to  form  RX2 ),  the  following  equation  is  obtained: 


F2  =/*+(/* 


"  f LOrx  )  y  LOrx  f SI 


(12) 


where  F2  is  the  frequency  of  the  signal  denoted  in  the  figure  as  RX2  and  is  given  with  respect  to  the 
system  timescale.  In  general,  the  values  fLOtx  and  fLOrx  are  tuned  such  that  their  difference  models  the 
nominal  value  of  fSLO  .  If  we  make  the  assumptions  that  the  satellite’s  LO  is  ideal  and  that  the  tuning  is 
perfect,  Equation  (12)  becomes 
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^2  ftx  y of  LOrx  * 


(13) 


Using  the  transformation  Equation  (5),  the  transmitted  frequency  fu  would  be  directly  measured  at  the 
receiving  end  (with  no  up/down  frequency  translation  or  satellite  relay)  as 
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Since  y0fwrx>yoftx >  observed  frequency  shift  (13)  is  greater  than  the  shift  expected  from  the 
transformation  (14)  by  a  factor  mnr ,  which  can  be  derived  as  follows: 
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(15) 


Note  that  the  right  side  of  Equation  (15a)  is  simply  Equation  (13)  transformed  through  the  use  of  the 
transformation  Equation  (5). 


Similar  to  the  derivation  of  Equation  (12),  proceeding  from  right  to  left  produces 

^1  ~  ftx  y ()ftx  ( f LOtx  J LOrx)  f  ofLOtx  f SLO  ’ 

which,  under  the  same  assumptions  as  before,  becomes 

U  ftx  y  0  (.ftx  f LOtx  ) 


(16) 


(17) 


and  it  is  seen  that  when  the  non-reference  station  is  the  transmitting  station,  the  fractional  frequency  bias 
is  multiplied  by  a  factor  of 


m=fwSL  +  l  (18) 

ftx 

where  mr  is  a  coefficient  of  the  clock’s  fractional  frequency,  as  seen  by  the  receiving  reference  station. 

From  the  above  discussion  we  see  that,  due  to  the  fractional  frequency  multiplication  effect  of  the  Carrier- 
Phase  Two-Way  system,  there  is  a  reduction  in  the  effective  measurement  quantization  and,  therefore,  a 
reduction  in  the  system’s  sensitivity  to  noise.  As  a  concrete  example,  consider  a  receiving  station  B  with 
a  clock  that  is  oscillating  1  mHz  slower  than  the  clock  at  a  transmitting  Station  A,  where  both  clocks  have 
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nominal  frequencies  of  5  MHz.  This  equates  to  a  fractional  frequency  deviation  of -2  x  10"10.  Suppose 
further  that  the  receiving  Station  B  is  capable  of  making  frequency  measurements  rounded  to  the  nearest 
tenth  of  a  hertz.  Now,  if  both  stations  are  collocated,  no  satellite  communication  is  necessary  and  Station 
B  can  take  a  direct  measurement  of  the  Station  A  output.  With  an  unlimited  resolution  and  in  the  absence 
of  measurement  noise,  Station  B  will  record  a  value  of  70,000,000.0140  Hz  (assuming  an  IF  nominal 
transmit  frequency  of  70  MHz).  However,  the  actual  value  recorded  is  rounded  to  70,000,000.0  Hz  and 
the  clocks  are,  therefore,  erroneously  determined  to  be  in  syntonization.  Now  consider  a  full  Ku-band 
satellite  transfer  (as  opposed  to  the  collocated  scenario)  where  the  transmitter  LO  and  the  receiver  LO  are 
tuned  to  model  a  satellite  turnaround  frequency  of  2  GHz,  with  the  nominal  values  of  13.93  GHz  and 
1 1.93  GHz,  respectively,  to  produce  a  14/12  GHz  transmit/receive  pair  when  mixed  with  the  70  MHz  IF. 
Using  the  coefficient  derived  in  Equation  (15),  the  fractional  frequency  of  the  received  signal  is 
multiplied  from  -2  x  10"10  to  approximately  -3.43  x  10"8  and  the  rounded  measurement  PB  is  then  equal  to 
70,000,002.4  Hz.  From  this  value,  the  two  Stations  will  deduce  that  their  fractional  frequency  deviation 
is 


_  L  ~Fb  _  (70,000,000- 70,000,002.4) Hz  ^  x 
70  mnrftx  171.43*  70,000,000  MHz 

and  the  result  is  very  close  to  the  actual  clock  frequency  deviation  (i.e.  -2  x  10"10). 


(19) 


IV.  SOLVING  THE  SYSTEM 


The  four  nonlinear  system  equations  (Equations  (6)-(9))  may  be  solved  by  several  different  numerical 
techniques.  One  method  begins  by  linearizing  the  four  equations  about  some  initial  estimation  vector 
X  =  [fsLoXX.df]  as  follows  [3]: 
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where 


SF,  =  F,  (f slo  ,ki,k2,df)-Fi  (fSLO  ,kl,k2,df). 


(21) 


For  simplicity,  Equation  (20)  is  re-written  as 


SF  =  GSx 


(22) 


where  G  is  generally  referred  to  as  the  geometry  matrix.  The  state  estimate  error-vector  Sx ,  of  the  above 
equation,  contains  updates  to  the  estimation  vector  x  and  summing  the  two  improves  the  accuracy  of  the 
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values  contained  in  5c .  This  method,  commonly  known  as  the  Newton-Raphson  method,  may  then  iterate 
while  Sx  converges  to  zero. 


However,  the  determinant  of  matrix  G  is  equal  to  0  and  the  linearized  system  is,  therefore,  singular.  As 
such,  no  unique  solution  exists  and  the  matrix  system  given  by  Equation  (22)  is  not  directly  solvable,  but 
necessitates  some  additional  numerical  method,  such  as  Singular  Value  Decomposition  (SVD).  However, 
it  is  known  that  the  values  kX9  k2,  and  df  will  be  very  close  to  0,  and  it  can  be  shown  that  the  system  is 

resistant  to  imperfect  approximations  of  fSLO  .  Therefore,  if  the  estimation  kx  =k2  =df  =  0  is  made,  and 
fSLO  is  assumed  to  be  known  and  constant  at  its  nominal  value,  a  single  iteration  will  suffice  and 
Equation  (20)  is  forced  to  the  static  equation  set 


^11  ~  ftx  +  f SLO 
^22  ~  ftx  +  fsLO 
F\2  ~ftx  +  fsLO 
f^2\  ~  ftx  +  fsLO 


"  f  SLO  +  ^ ^ 

0 

ftx 

-  f SLO  +  ftx 


0 

_  f SLO  +  2  ftx 
~  f SLO  ftx 

ftx 


0 

fsLO 

f 

J  sys 

f SLO  ~  ftx 

f 

J  sys 

fu 

f 

J  sys 


kx 
k  2 
df 


(23) 


If  Equation  (23)  is  further  reduced  by  treating  kx  as  an  independent  variable,  the  system  may  be 
generalized  by  the  geometric  figure  illustrated  below. 


Figure  5.  The  three  solution  points  and  their  average. 


As  shown  in  Figure  5,  Equation  (23)  produces  three  lines  with  three  different  points  of  intersection: 


df PI  =  f sys 
df  P2  ~  f  sys 
dfpi>  ~  f sys 
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V 
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(ftx  fsLO  )(^1 1  ^22  )  +  (f SLO  ^  ftx  )(Tj  1  Fn  ) 

-Iftx(ftx-fsLo) 


(2 fx  ~fsLo)(Fn -F2l)-ftx(Fn  ~F2S 

-  3  ftx  (ftx  -  fsLO  )  y 

f«(F2i-Fn)  +  fsLo(Fn-F2i)S 

Iftx(ftx-fsLo)  J 


\ 


(24) 


where,  if  the  equality  of  Equation  (10)  holds,  each  of  these  three  solution  points  is  expected  to  be  equal. 
This  allows  the  user  the  option  to  use  just  one  equation,  involving  any  three  of  the  Fmn,  thus  sparing  a 
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receiver  channel  for  unrelated  business.  Alternatively,  to  minimize  the  effect  of  noise,  the  three  points 
may  be  averaged  as 


df  A  VG  ~ 


f 

J  sys 


' fsLoiFn  -^.)-3  fa{Fn  -  F2i)' 

ftxiftx  ~  f SLO  ) 


(25) 


If  all  four  channels  are  used,  an  integrity  check  can  be  made  based  on  the  consistency  of  the  three 
solutions  in  Equations  (24). 


V.  RESIDUAL  ERRORS 

Our  previous  assumption  that  fSLO  is  static  and  nominal  will  cause  our  resultant  approximate  solution  to 
deviate  from  its  ideal  position.  Likewise,  inaccuracies  of  the  linearizations  (i.e.  truncation  of  the  high- 
order  terms  from  the  Taylor  Series  expansions)  will  cause  errors  in  the  solution.  Other  sources  of  error 
include  (but  are  not  limited  to)  troposphere  dynamics,  ionosphere  dynamics,  system  noise,  measurement 
noise,  and  multipath.  An  attempt  to  approximate  the  effects  of  some  of  these  error  sources  is  provided 
below  in  Table  1  [4,5]. 


Source 

Magnitude 

Resulting  Error  (ff)* 

Modem  quantization 

XmHz  *1.43x10-" 
It)  MHz 

~2xl0"13  /  tau 

Linearization 

k  =  6.7xl0"9 
fsLo  =  1x10 

dfciocks=  lxl  0"11 

<  lxlO16 

Unmodeled  zenith 
troposphere 

500  ps/day 

~  2.6xl0'15 

IGS  troposphere  maps 

10  ps/hr  RMS 

~  1.3xlO'15 

Unmodeled  zenith 
ionosphere  ** 

50  TECU  /  12hrs 

-  5xl0"15 

IGS  ionosphere  maps  ** 

3  TECU  /  2  hrs 

~  3xl0"15 

Locally  measured 
sinusoid  model  ** 

2  TECU  /  24  hrs 

<  lxlO"16 

*  Computed  for  Ku-band  uplink  and  downlink  frequencies  of  12  and  14  GHz 
**  At  GPS  frequencies,  2  TECU  causes  a  delay  of  roughly  1  ns 

Table  1.  Sensitivity  of  derived  frequency  transfer  to  modeling  errors. 


These  error  sources  are  non-Gaussian.  The  magnitudes  used  in  the  table  depict  “worst-case,”  but 
observable,  variations.  For  tau  <100  seconds,  the  dominating  error  is  the  1  mHz  quantization  of  the 
TWSTFT  modems,  which  currently  provides  1  mHz  resolution  at  IF.  Future  modem  revisions  are 
expected  to  improve  this  figure. 


36th  Annual  Precise  Time  and  Time  Interval  (PTTI)  Meeting 


VI.  RESULTS 

Two  preliminary  tests  have  been  prepared  for  inclusion  in  this  report;  the  results  of  these  test  runs  make  it 
clear  that,  while  the  theoretical  aspects  of  the  system  are  sound,  more  work  is  required  before  carrier- 
phase  TWSTFT  can  become  fully  operational. 

The  first  test  case  was  conducted  using  two  collocated  stations  at  USNO.  The  86,400  fractional 
frequency  data  points,  shown  in  Figure  6  (below),  was  calculated  with  the  dfPZ  result  of  Equation  (24). 
The  data  show  a  modulating  frequency  with  a  rate  of  approximately  0.8  mHz,  or  about  21  minutes  (see 
Figures  7  and  8),  and  this  is  reflected  in  the  Allan  deviation  plot  shown  in  Figure  9.  Subsequent  tests 
have  suggested  that  the  problem  may  be  traced  to  low-level  frequency  generations  within  the  spread- 
spectrum  modems.  Specifically,  if  a  frequency  which  is  to  be  generated  within  a  modem’s  direct  digital 
synthesis  (DDS)  block  is  not  an  exact  multiple  of  the  DDS’s  minimum  incremental  resolution,  it  is 
compensated  for  by  the  phase-lock  loop  and  an  artificial  frequency  modulation  can  occur  at  the  PLL 
output.  However,  it  is  also  possible  that  the  witnessed  data  variations  are  simply  due  to  temperature 
fluctuations  or  other  environmental  effects.  These  possibilities  are  currently  under  investigation. 

Figure  10  shows  the  Allan  deviation  from  a  day  of  carrier-phase  TWFTFT  sessions  involving  stations  at 
both  USNO  and  Physikalisch  Technische  Bundesanstalt  (PTB).  While  we  were  unable  to  produce 
enough  uninterrupted  data  to  provide  stability  calculations  beyond  4000  seconds,  it  is  evident  that  there 
exist  oscillations  within  these  data  as  well.  It  is  interesting  to  note,  that  the  oscillations  in  this  data  set  are 
at  a  higher  rate  than  the  0.8  mHz  frequency  seen  in  the  collocated  data. 


VII.  CONCLUSION 

In  this  paper,  we  have  refined  previously  published  techniques  for  calculating  the  carrier-phase  two-way 
solution.  The  paper  suggests  a  method  of  solving  the  system  and  then  presents  a  simple  error  analysis. 
Worst-case  error  dependencies  related  to  quantization,  linearization,  and  atmospheric  effects  of 
ionosphere  and  troposphere  are  shown.  We  conclude  that,  with  an  ideal  carrier-phase  modem,  a  2  x  10"15 
frequency  transfer  at  1  second  is  possible. 

Two-Way  modems  used  at  many  of  today’s  TAI  laboratories  have  limitations  due  to  the  modems’  current 
measurement  resolution.  Our  bench  tests  and  live  experiments  demonstrate  these  limitations.  While  it  is 
unclear  what  is  causing  the  additional  noise,  several  possibilities,  including  frequency  synthesis  within  the 
modems,  are  being  explored.  These  unknown  effects  must  be  well  understood  before  the  carrier-phase 
TWSTFT  system  becomes  fully  operational,  in  order  to  maintain  consistency  and  accuracy  within  the 
system’s  error  budget. 
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Figure  6.  Fractional  frequency  collocated  data  (1-second  data  points)  collected  at  USNO.  Outliers  that 
exceed  a  sigma  factor  of  10.0  were  removed. 
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Figure  7.  A  40-point  running  average  of  the  collocated  fractional  frequency  data  accentuates  the 
modulating  0.8  mHz  frequency. 
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Figure  8.  The  power  spectrum  of  the  collocated  fractional  frequency  data  shows  that  the  frequency  of  the 
modulating  signal  is  approximately  0.8  mHz. 


Dote;  01/11/05  Time;  11:52:4  6 


Doto  Point 4  1  thru  86400  of  342821 


Tou  1.0000000e400  Hie:  53138. output  _vio3ChonnelSoluti 


FREQUENCY  STABILITY 

Collocated  at  USNO  (86400  data  pnts) 


Tqu 

Sigmo 

!  I.OOe+OO 

4.14e-13 

.  2.00e+00 

2.35e-13 

_ 

4.00e+00 

I.23e-13 

8.00e  +  00 

6.8le-l4 

•  1.60+01 

3.85e-14 

- 

3.20+01 

2.60-14 

6.4-Oe  +  Ol 

3.l6e-l4 

1.28e+02 

5.07e-14 

: 

2.56e+02 

8.02e-14 

.j 

5.12e+02 

1.00e-13 

- 

1.02e+03 

3.67e-14 

- 

2.05e+03 

2.94<?-14 

4.i0e+03 

i.48e  14 

f 

,  8.i9e+03 

I.l8e-14 

,  1.64e+04 

8  28e-15 

101  10z  103  104 

Averaging  Time,  x,  Seconds 


Figure  9.  Alan  Deviation  plot  of  the  collocated  results  (1 -second  data). 
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Figure  10.  Carrier-Phase  TWSTFT  session  between  USNO  and  PTB. 
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Questions  and  Answers 


TOM  BAHDER  (U.S.  Army  Research  Laboratory):  When  you  set  up  the  whole  model,  what 
reference  frame  -  are  you  in  a  rotating,  or  are  you  in  an  ECI  frame? 

BLAIR  FONVILLE:  Well,  the  satellite  is  geostationary,  so  we  set  up  each  equation  and  derive  all 
frequencies  with  respect  to  a  system  reference  station. 

BAHDER:  But  then  you  are  implementing  everything  in  a  rotating  frame,  right?  And  my  other  question 
is  then  there  must  be  some  kind  of  a  Sagnac  effect  that  you  should  be  seeing. 

FONVILLE:  Yes,  we  may  see  some  slight  drifts  due  to  uncorrected  Sagnac  in  the  long  term,  but  we 
haven’t  been  able  to  take  data  over  15-minute  spans.  Like  I  say,  diurnal  effects  are  also  something  we 
need  to  investigate. 

BAHDER:  You  could  avoid  the  Sagnac  term  if  you  did  the  whole  implementation  in  the  ECI  frame.  But 
then  you  would  have  relative  velocities  of  the  ground  stations  and  the  satellites,  and  that  would  have  to 
come  into  the  problem. 

FONVILLE:  The  relative  velocities  of  the  ground  stations  to  the  satellite  is  something  we  do  account 
for.  Thank  you  for  that  comment. 

CHRISTINE  HACKMAN  (University  of  Colorado):  You  put  up  a  viewgraph  where  you  were  talking 
about  all  the  errors,  and  you  showed  two  parts  in  1015  for  the  troposphere. 

FONVILLE:  Yes,  I  knew  you  were  going  to  ask  something  about  that  -  after  your  presentation. 

HACKMAN:  So  are  you  assuming  it  is  not  going  to  cancel?  Could  you  elaborate  on  that  a  little? 

FONVILLE:  The  troposphere  will  not  cancel  because  we  are  operating  at  two  different  sites. 

HACKMAN:  I  guess  what  I  normally  understood  in  regular  Two-Way,  you  are  kind  of  going  for  this. 

FONVILLE:  Oh,  I  see.  Yes,  well,  we  have  a  transfer  from  one  site  to  the  other  and  vice  versa,  but  we 
also  have  two  loopbacks.  And  so  those  two  loopbacks  are  where  it  does  not  cancel.  We  have  a  loopback 
on  one  side,  we  have  a  loopback  on  the  other;  they  have  different  tropospheric  paths,  so  it  is  not  going  to 
cancel  there. 

WLODZIMIERZ  LEWANDOWSKI  (Bureau  International  des  Poids  et  Mesures):  I  have  a  short 
question.  What  is  your  feeling  about  the  future  of  this  approach?  When  can  we  expect  the  operation 
technique? 

FONVILLE:  That  is  a  good  question.  This  is  all  very  preliminary,  but  the  results  are  encouraging.  So  it 
does  give  us  some  motivation  to  continue  the  effort.  It  can  be  in  operation  relatively  quickly.  But  it  is 
contingent  on  how  the  longer  periods  look.  For  instance,  we  have  had  some  data  that  -  well,  in  the  15 
minutes  that  I  showed,  it  seems  to  be  dropping  without  any  problem.  But  we  have  concatenated 
contiguous  sets  of  data,  and  tried  to  do  some  numerical  calculations  on  it.  And  it  seems  to  level  off. 
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So  we  may  still  have  some  processes  that  we  have  not  yet  fully  investigated.  So  that  is  a  difficult 
question  to  answer. 


